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Wideband Reference-Plane Invariant Method for 
Measuring Electromagnetic Parameters of Materials 

Khattiya Chalapat, Kari Sarvala, Jian Li and Gheorghe Sorin Paraoanu 



Abstract — This paper presents a simple and effective wideband 
method for the determination of material properties, such as the 
complex index of refraction and the complex permittivity and 
permeability. The method is explicit (non-iterative) and reference- 
plane invariant: it uses a certain combination of scattering 
parameters in conjunction with group-velocity data. This tech- 
nique can be used to characterize both dielectric and magnetic 
materials. The proposed method is verified experimentally within 
a frequency range between 2 to 18 GHz on polytetrafluoroethy- 
lene and polyvinylchloride samples. A comprehensive error and 
stability analysis reveals that, similar to other methods based on 
transmission/reflection measurement, the uncertainties are larger 
at low frequencies and at the Fabry-Perot resonances. 

Index Terms — Reference-plane invariant method, refractive 
index measurement, electric permittivity, magnetic permeability, 
transmission/reflection method, broadband measurement. 



I. Introduction 

IN recent times, fast and accurate knowledge of the electro- 
magnetic properties of materials in the microwave range is 
increasingly required in the design and development process 
of a vast number of industries, spanning from food processing 
to communication systems. Also in basic science research, 
understanding and measuring material parameters, such as the 
complex refractive index n and the complex permittivity e, is 
an important fundamental task. 

Nonresonant techniques such as transmission and reflection 
measurements are largely used nowadays for characterizing the 
electromagnetic properties of materials [1]; the fundamentals 
of these techniques have been already established in the 1970's 
by the seminal papers of Nicolson, Ross [2] and Weir Q. 
These techniques are relatively simple and accurate. They 
have the advantage of broadband characterization of materials 
and devices. For almost four decades they have been widely 
applied to measure the permittivity and permeability of various 
synthetic and natural materials. 

A known drawback of the original Nicolson-Ross-Weir 
(NRW) method is that it requires the transformation of S- 
parameter measurements from the calibration reference planes 
to the surfaces of the material. The phases of the transmission 
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and reflection signals are strongly dependent on the positions 
of the reference planes, so the uncertainties in the transforma- 
tion of S-parameters can add significant errors. The precision 
of this transformation can be enhanced in various ways, for 
examples by adding more steps to the calibration process 
or running extra calculation algorithms which complicate 
the measurement. Besides the Nicolson-Ross-Weir algorithm, 
other methods based on transmission/reflection measurements 
exhibit the same sort of difficulties. 

An important step forward has been achieved in 1990, when 
Baker-Jarvis and collaborators showed that it is possible to de- 
rive S-parameter equations which are reference-plane invariant 
|4]. Using some of these equations, they showed that it is pos- 
sible to extract the values of the material parameters by using 
an iterative algorithm. This algorithm requires as input some 
initial values for permittivity and permeability, and therein lies 
one of its limitations: a good guess is needed, otherwise the 
algorithm can produce wrong results. Although valuable work 
has been done recently to improve these ideas |5)- iflOl . there is 
currently no universally-accepted best method. Each proposed 
technique has a number of advantages and disadvantages and 
the weight to be attached to each depends on the specific 
application. 

The purpose of this paper is to present a methodology 
for extracting the complex material parameters from trans- 
mission/reflection measurements which combines ideas from 
the NRW and the Baker-Jarvis techniques. More precisely, 
the scattering parameters are combined into a specific set of 
reference-plane invariant equations (similar to Baker-Jarvis), 
and the equations are used together with group velocity data 
(similar to NRW) to obtain the complex permittivity and 
permeability. Surprisingly, this results in a simple, explicit, 
and reference-plane invariant methodology which can be used 
to characterize both dielectric and magnetic materials. With 
respect to NRW, the advantage of our method is that it uses 
reference-plane invariant quantities, therefore the errors due 
to calibration to the two material-air interface are eliminated. 
With respect to the Baker-Jarvis algorithm, the improvement 
consists in the use of additional information about the sample, 
extracted from group velocity measurements. This removes the 
ambiguity in the determination of the phase. Also, since our 
results do not depend on choosing good initial values for e 
and /i as in Baker-Jarvifl dielectric materials with unknown 
properties and also materials with magnetic properties at high 
frequency can be characterized. 

'The method of Baker-Jarvis is tailored for dielectric materials, thus one 
takes fjt r = 1 as initial value. 
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II. Theory 

In this section, the measurement is modeled within the 
framework of classical electrodynamics. We present a new 
algorithm which is reference-plane invariant and show how 
it can be used to determine the complex refractive index and 
the complex permittivity and permeability. 

We start by deriving the mathematical relations between 
the S -parameters and the material parameters. We describe 
the scattering of electromagnetic waves based on the multiple 
reflection model shown graphically in Fig. [TJ Within this 
model, the total reflection and transmission coefficients can be 
calculated using the superposition principle. Since a transverse 
electromagnetic (TEM) wave traveling a distance L picks up 
a phase change of 2irL/\, where A is the wavelength in that 
region, the propagation factor of the TEM wave traveling 
through the material of length L, as shown in Fig. Q]is given 

by 

P = e" 72i , (1) 
where 72 = iuj/v 2 = iujn 2 /c = iuj^/JJ^e^. 




Fig. 1. The model of multiple reflection between two interfaces. 

The total reflection coefficient is 

r tot = r + t 21 t 12 tp 2 + r 21 T 12 r 3 p 4 + . 

r(i - p 2 ) 
1 - r 2 P 2 ' 

1 - 
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where 
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Similarly, the total transmission coefficient in terms of V 
and P is 

_ p(i-r 2 ) 

Ptot — j - p 2 p 2 • ^ 

The standard model of a TR measurement is described in 
Fig. |2] The transmission line is divided into three regions. 
Typically, the regions of lengths Li and L 2 are assumed to 
be filled with air and the middle region of length L is filled 



with a material of relative permittivity e r = e/e and relative 
permeability \i r = fi/fio- The complex refractive index of the 
material is n = ^//.i r e r . We then take the permittivity, ej., and 
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Fig. 2. The model of a transmission line containing a material of length 
L. Lj (j = 1,2) represents the distance from the reference plane of the 
S-parameter measurement to the corresponding interface between air and the 
material under test. 

the permeability, /ii, to be equal to the permittivity of free 
space, £q, and, respectively, the permeability of free space, 
fio, and write Eq. (O in the form 



(6) 



z + 1 



where z — y/ 'n r /e r is the impedance relative to vacuum (the 
total impedance Z = J fi/e = Zqz, where Zq = J Jiq / eo = 
1207rf2 is the vacuum impedance). 

The determination of these two quantities, z and n, from 
the experimental data will be the main focus of the remaining 
part of this section. Based on the multiple reflection model, 
the S-parameters are expressed in terms of V and P as follows 

P 2 ) 



r tQt = e - 2 ^^ 



-271L 



2 r t , 



1 - r 2 P 2 : 
_ 2 ^ r(i-p 2 ) 

1 _ p2p2 ' 



(7) 



(8) 



and 



021 — 012 — e 1 tot — e — j, 2 p2 ■ 

(9) 

where 71 = iujni/c w iui/c (c is the speed of light in the 
vacuum). 

When there is no sample inside the transmission line, T = 0, 
72 = 71 and therefore 

S° t = e -7i(ii+i2+i) = e -7i(ia ir )_ (10 ) 

Eq. ( [Tot allows us to experimentally determine the airline 
length, L a j r , by calibrating the vector network analyzer and 
then measuring the transmission through the empty air line to 
obtain the phase of S 21 . 

In the next step, F and P are expressed in terms of 
the S-parameters. This is similar to the Nicolson-Ross-Weir 
algorithm, with the essential difference that neither T nor 
P depends on L\ and L 2 . Indeed, for airlines operating at 
relatively high frequencies, measurements of L\ and L 2 are 
prone to relatively large uncertainties. These errors will further 
propagate in the phase factors of the S-parameters, 



6 (2 7l Lj 



3 €{1,2}, 



(11) 
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leading to the larger errors of the phase factors at higher 
frequencies. 

Fig. [3] shows the schematic model of a reference-plane 
invariant measurement. We start by defining two quantities, 
which are related to measurable quantities, namely 

r 2 (i-p 2 ) 2 



A = 



22 



SuS: 

(1-r 2 ) 2 W 



(12) 



and 



B = e 2 ^ L -- L \S 21 S 12 - S n S 22 ) = Y~f^2- (13) 



Calibration Plane 



Calibration Plane 





£ 2 ^2 




L 



Fig. 3. Reference-plane invariant measurement model. In the diagram, L 
represents the length of the sample and L a ; r represents the distance between 
the calibration planes of the S-parameter measurement. 

Experimentally the 5-parameters are measured by a vector 
network analyzer (VNA), the airline length, L a h, is found via 
Eq. ( [Tol l, and the length of the sample, L, is measured before 
inserting the sample into the airline. 

Solving Eq. CE3 for P 2 , 

r 2 



P 2 = 



D 



1 + BY 2 ' 

and substituting back into ( fT2b . we obtain 



A = 



T 2 {\-Bf 



(b + t 2 )(i + bt 2 Y 

Eq. (Q3) can be solved to find T 2 

2 = -A(l + B 2 ) + (1-B) 2 
2AB 



(14) 



(15) 



(16) 



± 



-AA 2 B 2 + [A{l + B 2 )-{\- B) 
2AB 



where the sign in this equation is chosen so that |T| < 1. 
These expressions for P 2 and T 2 are manifestly reference- 
plane invariant. 

A very useful, simpler expression for P can be obtained if 
we define another quantity, R, directly related to the scattering 
parameters, 

S21 _ e< lL P(l - r 2 ) 
si\ ~ 1 - P 2 r 2 ' 



P = 



(17) 



Then Eq. ( fTTl ) can be solved for P. Substituting P 2 from < fT~4T > 
into the denominator of ( fTTl ), we obtain 



P = R 



1 



1 + pr 2 



-liL 



(18) 



From ([]]), assuming free space on either side of the sample, 
the complex refractive index can be determined by 



7i L 



■In 



or, more explicit, 



1 



71 L 



■In 



1 



1 + pr 2 



R 



(19) 



(20) 



The logarithmic function in ( TT9l > and d20l > is a multi-valued 
function, which results in an infinite number of discrete values 
for n . The physically correct solution must be chosen from 
these values. One way to do so is to check whether a chosen 
solution gives correct values for another measurable quantity 
or not. Also, this measurable quantity should not depend on 
L\ and L 2 . In the following, we will show that an appropriate 
such quantity is the group delay in the line. 

By definition, the group delay is a measure of a pulse signal 
transit time through a transmission line. The transit time of a 
wave packet is defined as 



x 
v„ 



(21) 



where x is the transit length and v g is the group velocity of 
the wave pulse. In this case, 



and 



L- 



df 



fn 



(22) 



(23) 



where r° is the group delay through an empty line, and T g is 
the group delay through the line with an inserted sample of 
length L. By comparing the calculated group delay with the 
measured group delay, 



.(measured) 



(calculated) 



0, 



(24) 



the correct refractive index can be determined. Another quan- 
tity which can be used as an alternative to the group delay is 
the group delay relative to the empty air-line, 



L ( An 



(25) 



derived from supstracting Eq. ( f23l from Eq. ( f22l . 

In the final part of this section, we describe how to extract 
the material parameters e r and fj, r . A situation of practical 
interest is the case in which the experimentalist already has 
some information about the material. For example, if chemical 
analysis provides additional proof that the material does not 
contain magnetic elements, one can take /i. r = 1 and determine 
the permittivity from e r = n 2 . As we will show later, this leads 
to better accuracies than the more general method presented 
below, which requires the determination of z. 

In many situations however, especially concerning materials 
under research which contain magnetic elements, magnetic 
impurities, or magnetic nanoparticles, it is not possible to 
know beforehand what the electromagnetic properties are. In 
these situations, one needs to use not only n but also the 
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relative impedance z; the material properties are then obtained 

as 

fx r = nz, e r = n/z, (26) 

where 



The reflection coefficient T is determined from Eq. ( TToT l. But 
this equation gives T only up to a sign, since ±r both satisfy 
Eq. (T5[ . To get the correct sign for T, we have to go back 
to Eq. (O or (O and check which one of ±T satisfies them. 
Note that this does not bring in additional errors, since it is 
just a sign check. 

It is useful to note that even without such a check, only 
minimal information about the properties of the material may 
be sufficient. Suppose that T is the correct solution leading to 
the correct set of material parameters e r and fi r . The properties 
of the conformal mapping, Eq. d27| >, imply that the opposite- 
sign solution — T corresponds to a relative impedance z _1 . The 
effect in the final result Eq. (|26| i is therefore simply to swap 
the values of permittivity and permeability. In many practical 
situations, an experienced experimentalist could recognize 
easily, given two complex numbers and minimal information 
about the chemical composition of the material, which one is 
the permittivity and which one is the permeability. 

III. Experiments 

To verify the proposed method, the complex refractive 
indices of polytetrafluoroethylene (PTFE) and polyvinylchlo- 
ride (PVC) were determined. The reference plane invariant 
algorithm was tested by measuring the PTFE and PVC sample 
at various positions relative to the calibration planes of the S- 
parameter measurements. The group delays and S-parameters 
were measured by using an Anritsu 37369D vector network 
analyzer. The measurement setup is shown diagramatically 
in Fig |4] Transverse electromagnetic waves are transmitted 
between the measurement ports via a 7mm precision air 
line set which consists of an outer conductor with beadless 
connectors and a 3mm center conductor. The air line set was 
connected to the VNA ports by Anritsu 34ASF50-2 female 
adapters. 

Vector Network Analyzer 
40.0 MHz - 40.0 GHz 

? X 

7 mm coaxial air line, DC to 1 8 GHz 




Calibration Plane Calibration Plane 



Fig. 4. The measurement setup used to measure the complex refractive 
indices of test samples. The 7-mm coaxial air line was used as the sample 
holder to conduct transverse electromagnetic waves. From the diagram, an 
inserted sample can be seen in the middle of the air line. 



The full 12-term calibration, excluding isolation, was per- 
formed prior to the measurements. First the empty air line 
was measured to obtain 5fi and t°, from which the total 
length of the air line L a ir between the calibration planes was 
inferred to be 17.3193 cm. Then a toroidal sample was inserted 
between the inner and the outer conductor of the air line, and 
the measurement was repeated again. 

Measurements on a 20.00mm PVC sample give the results 
shown in Fig. [5] and [6] Within the whole range of measured 
frequencies (40.0 MHz to 18.0 GHz) there are 450 sampling 
points. Most of the samplings gave results which are close to 
the average values. For example, at a frequency of about 12 
GHz, the measurement gave n = 1.595 - 0.012 i. Comparing 
the real values, n', and the imaginary values, n", of the 
measured refractive index, the algorithm leads to very stable 
results for the imaginary part, but the real part contains two 
discontinuities at 4.60 GHz and 14.04 GHz. These frequencies 
correspond to integer multiples of one-half wavelength inside 
the sample. By ignoring the discontinuity, the real part of the 
refractive index can be interpolated as 1.605 around 4.60 GHz. 
For a wavelength of twice the sample length, this corresponds 
to a frequency of 4.67 GHz, close to the measured value of 
4.60 GHz. 
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Fig. 5. Multiple roots of Eq. U9t which were used to determine the real 
part of the refractive index shown in Fig. [6] 

Since PVC is non-magnetic, the relative magnetic per- 
meability is approximately equal to 1. We can therefore 
use the measured refractive index to determine the electric 
permittivity, e r = n 2 . Fig. [7] shows the measured result for 
the PVC sample of length 20.00 mm obtained by using the 
proposed method and the Nicholson-Ross-Weir method. We 
can see that for a non-magnetic material, such as PVC, the 
complex permittivity can be measured using this method with 
relatively high accuracy. 

To further examine and verify the proposed algorithm, we 
also performed a measurement on a PTFE sample of length 
20.00 mm. The results are shown in Fig. [8] At 10 GHz, the 
index of refraction has the value of about 1.415 - 0.004i. The 
spectrum of the real values, n', is discontinuous at 5.24 and 
15.84 GHz, while the spectrum of the imaginary values, n", 
is relatively stable over the whole range of frequency, but we 
can see one relatively large peak near 15.84 GHz. In Fig. [9] 
the electric permittivity of PTFE is determined from e r = n 2 . 
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Fig. 6. The complex refractive index of a 20.00 mm PVC sample determined Fig. 8. The complex refractive index of a 20.00 mm PTFE sample determined 
by using the reference-plane invariant method. by using the reference-plane invariant method. 
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Fig. 7. Comparison of the complex permittivity obtained using the Nicholson- 
Ross-Weir method (NRW) and our reference-plane invariant method (RPI) for 
a PVC sample of length 20.00 mm. The reference-plane positions used in the 
NRW algorithm are measured by the vector network analyzer. 



The spectra obtained by the new method show much less errors 
compared to the ones from the NRW algorithm, especially at 
higher frequencies where we have larger uncertainties in L\ 
and L 2 - 

The method is invariant with respect to the position of 
the reference-planes, therefore, theoretically, the position of 
the sample inside the air line should not affect the results. 
To verify this statement, a number of measurements were 
repeated by placing the samples at different positions inside 
the air line. The statistical analysis of the measured results are 
shown graphically in Fig. [iOlflTI The error bars represent the 
uncertainties, i.e. standard deviations, of the measurements. It 
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Fig. 9. Comparison of the complex permittivity obtained using the Nicholson- 
Ross-Weir method (NRW) and our reference-plane invariant method (RPI) for 
a PTFE sample of length 20.00 mm. The reference-plane positions used in 
the NRW algorithm are measured by the vector network analyzer. 



can be seen that the uncertainties of either n' or n" are lower 
at higher frequencies. 

For PTFE, the uncertainties of the real values vary slowly 
with frequency, while the uncertainties of the imaginary values 
are lower at higher frequencies and become larger again near 
18 GHz. This might be explained in terms of energy losses 
within the line because the measurement on the empty line 
shows higher losses at frequencies near 18 GHz. In PVC, this 
effect is not observed because the material losses are larger 
and tend to dominate the air line imperfections. 



7 



1.63 



1.62 



K 1.61 



1.60 



1.59 



PVC 20.00 mm 



2 4 6 8 10 12 14 16 18 
/ [GHz] 



0.00 



-0.02 



-0.04 



-0.06 



PVC 20.00 mm 



8 10 12 
/ [GHz] 



14 16 n 



Fig. 10. The refractive index of PVC measured by using a sample of length 
20.00 mm. 



IV. Error Analysis 

Using the proposed method, the complex refractive index is 
determined by measuring the transmission/reflection signals, 
the group delay, the air line length and the sample length. 
According to the rules of error propagation, higher errors 
in the sample length and S-parameter measurements will 
lead to higher uncertainties in the measured refractive index. 
Experimentally, the transmission and reflection signals are 
measured by a VNA, therefore the accuracy of the measure- 
ment is limited somehow by the VNA uncertainties. In fact, 
there are also other causes of errors to consider, such as the 
imperfections of the air line and the air gaps between the 
surfaces of the sample and the line. The errors due to the air 
gaps can be compensated by mathematical models. Equations 
for the air gap correction can be found in the literature liTTTl . 

ma. 

The analysis is simplified by assuming that the errors due 
to air gaps, connector mismatches and air line losses are 
small compared to the errors due to the sample length and S- 
parameter measurements. The uncertainties of reference-plane 
positions are also neglected because the final equations are 
theoretically reference-plane invariant. 

Experimentally, the S-parameters are independently mea- 
sured by a VNA, so if the errors of the magnitudes and 
phases of the S-parameters are assumed to be independent, 
the relative uncertainty can be calculated using the rules of 
error propagation as follows: 
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The index a € {11, 12, 21, 22, 21o} runs over the scattering 
parameters, S a — \S a \exp(id a ). (The index 21o refers to 



the transmission through the empty air line.) Eq. 
further expressed as 
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Fig. 1 1 . The refractive index of PTFE measured by using a sample of length 
20.00 mm. 



In the remaining part of this section, the errors involved 
in our experimental determination of the dielectric properties 
of PVC and PTFE will be analyzed. The derivatives in d29l ) 
can be analytically calculated, but the results are complicated 
and it is in fact simpler to numerically calculate the relative 
uncertainties in (l28l i. Fig. fT2land Fig. [T3l shows the numerical 
results for two models, one with n = 1.595 - 0.012i and another 
with n = 1.416 - 0.003i. 

The numerical calculations are done based on the assump- 
tion that the materials are non-magnetic. That is, the relative 
permeability of each model is taken to be 1, fi r = 1, and 
the relative permittivity is equal to the square of the refractive 



index, e r — 



The errors due to the sample length and S- 



parameter measurements are set as follows: SL — 0.1 mm, 
S\S a \ = 0.002, 56 n = S9 22 = 3° and 66 21 = S9 12 = S9 21o = 
1°. The error SL is set greater than the actual value in order 
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Fig. 12. The relative uncertainty of refractive index for a model sample of 
length 20.00 mm and n = 1.595 - 0.012i. The middle graph shows the same 
result as the upper graph at a magnified scale. The relative uncertainties of 
the real part, 8n' /n', are very high around the resonant frequencies, but close 
to zero elsewhere. Both the real and imaginary spectra show relatively low 
uncertainties at high frequencies. 

to investigate the behavior of the corresponding uncertainty. 
In the experiment, the PVC and PTFE samples were precisely 
machined, and the uncertainties of the lengths and diameters 
of the samples are less than 0.04 mm. 

It can be seen from the relative uncertainties in Fig. [12] 
and the experimental results in Fig. [6] and Fig. [10] that the 
errors determined theoretically correspond well with the errors 
determined experimentally. As in other methods based on 
transmission/reflection measurements, uncertainties are lower 
at higher frequency. This is understandable, because a signif- 
icant change in the phase due to the presence of the material 
in the air line requires that the sample length should not be 
too small compared to the wavelength. Around 4.70, 9.39 



Fig. 13. The relative uncertainty of refractive index for a model sample of 
length 20.00 mm and n = 1.416 - 0.003i. The middle graph shows the same 
result as the upper graph at a magnified scale. The imaginary spectrum shows 
some small peaks at resonant frequencies. 

and 14.08 GHz, the simulation shows discontinuities in the 
real part spectrum. As mentioned before, the origin of these 
discontinuities is of physical nature, namely the occurrence of 
Fabry-Perot resonances in the sample. However, our results 
for n (and also for e r — n 2 in the case of non-magnetic 
materials) are less sensitive to this effect than those obtained 
via the original NRW algorithm: this can be seen already in 
Fig. [7] The difference is mathematical: in the NRW algo- 
rithm the discontinuity is generated by the relative impedance 
z = (1 — r)/(l + T) |6], while for non-magnetic materials our 
algorithm avoids calculating this quantity. A comparison with 
the situation in which z needs to be calculated is given also 
below in Section [V] 

Fig. Q~3] shows the simulated results for a sample with 
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Fig. 14. The relative uncertainty of n" for a model sample of length 
20.00 mm and n = 1.416 - 0.003i. Each line represents relative uncertainty, 
(1/n") X (8n" /&x) X Sx, caused by a single-source error, 8x. The upper 
graph shows the relative uncertainty due to the S-parameter errors, i.e. 
x 6 {611,621, |>Sii|, IS12I, 1^21 1}- The lower graph shows the uncertainty 
due to the sample length error. 



V. Theoretical considerations on measurement 

ERRORS FOR VARIOUS TYPES OF MATERIALS 

For dielectric materials, we have already analyzed measured 
data for PVC and PTFE; for magnetic materials, no reliable 
standards or reference materials have emerged. In this section, 
we will give a systematic analysis of measurement errors for 
both dielectric and magnetic materials based on simulated 
experiments. First, a set of S -parameters and group delays 
are generated based on given values of the sample length L, 
the air-line length L a j r , the complex electric permittivity e r 
and the complex magnetic permeability fi r . Then errors are 
added to these data and the results are used as the inputs of 
the reference-plane-invariant algorithm. 

A. Low-loss dielectric materials 

Figs. [l5lfT6l show the results of simulated measurements of 
the refractive index for two zero-loss dielectric materials (e r 
= 5 and 7) with sources of errors S\S a \ = —0.001, S9 a = 
1° and 8L a i r = 0.1 mm. The length of the sample was taken 
L = 5 mm. 
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n = 1.416 — 0.003i. It is clear from the simulation that 
the relative uncertainties in n' and n" are again prone to be 
higher at lower frequencies. The uncertainty in n" oscillates in 
frequency domain and has small peaks at resonant frequencies, 
i.e. frequencies corresponding to integer multiples of one-half 
wavelength inside the sample. 

Fig. Q3] shows a comparison between the uncertainties due 
to each source of error. The uncertainties originating from the 
S-parameter errors oscillate and decrease in amplitude when 
the frequency is increased. This explains why the total relative 
uncertainties of n" shown in Fig. Qj] behave in a similar way. 
The error in the sample length measurement leads to high 
uncertainties at frequencies corresponding to integer multiples 
of one-half wavelength, but very low uncertainties at other 
frequencies. In experiments, sample lengths can be measured 
with high accuracy, so the corresponding errors are much 
smaller than the simulated results. 

We have also simulated the relative uncertainties of n' and 
n" caused by errors in the air-line length, and found that the 
corresponding uncertainties behave similarly to those caused 
by the sample-length error. Experimentally, the error due to the 
air-line length depends on the accuracy of the TRU calibration. 



Fig. 15. Simulated real-part spectra of the refractive indices for dielectric 
materials with e r = 5 and 7. 
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Fig. 16. Simulated imaginary-part spectra of the refractive indices for 
dielectric materials with e r = 5 and 7. 



Similarly to what we have shown in Fig. [14] we see 
clearly that when the sample length is small compared to the 
wavelength, the errors become significantly higher at lower 
frequencies. We also check the code by setting S\S a \ = 0, 
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56 a = 0°, 6L air 
initial values of e r 



7.6 



0, and find that we recover the exact the NRW algorithm: measurement errors become magnified 

at and around the frequencies at which the reflection goes to 
zero. The difference is that here we do not have errors due to 
the reference-plane positions. 



7.4- 



7.2 



7.0 



6.8 





L = 5 mm 


\ 8|SJ, 89 






5L . 




/ " 






/ 

Exact (input) 





2 4 6 8 10 12 14 16 IS 
/ [GHz] 



Fig. 17. Simulated real spectra of electric permittivity for a dielectric material 
with e r = 7, calculated from the relation: e r = n 2 . The errors are set as 
follows: S\S a \ = -0.001, 58 a = -1° and SL air = 0.1 mm. 
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Fig. 18. Simulated real spectra of electric permittivity and magnetic 
permeability (extracted as e r = n/z and fi r = nz) for a dielectric material 
with e r = 7, with the errors in the S-parameters (upper plot) and L a i r (lower 
plot). 

In Fig. [T7] we present the permittivity spectra determined 
from the complex refractive index using e r — n 2 , and in 
Fig. [18] we show the same spectra obtained using the version 
of the algorithm in which e r and fi r are extracted via the 
relative impedance z, e r — nj z and /j, r = nz. In the second 
case, the errors around the resonant frequencies get larger. 
The mathematical origin of these errors is similar to that of 



B. Lossy dielectric materials 

For lossy dielectric materials, resonances do not lead to zero 
reflected signals at the resonant frequencies, which results in 
smaller errors in the permittivity and permeability spectra. 
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Fig. 19. The simulated spectra of the complex refractive indices for a lossy 
dielectric material with e r = 5 — 0.2i, /j, r = 1. The errors are set as follows: 



S\S a \ 
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-1° and 8L„ 
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To quantify how the energy loss affects the errors, the S- 
parameters generated for e r — 5 — 0.2i are used as the input 
of the reference-plane invariant algorithm. The simulation is 
done with <5|5 a | = -0.001, 59 a = -1° and 5L air = 0.1 mm. 
The results are shown in Fig. [T9j indicating that the errors are 
quite low over the entire frequency range. It is now possible 
to extract, with small errors, e r by using e r = n 2 . 

Alternatively, one can extract the complex e r using the 
relative impedance z, with results shown in Fig. |20] Similar 
to zero-loss materials, measurement uncertainties cause sig- 
nificant errors at and around Fabry-Perot resonant frequency, 
i.e. frequency corresponding to integer-multiple of one-half 
wavelength inside the sample. However, for lossy materials, 
only finite peaks are found, with no divergence, in both spectra 
of e r and fi r . 

C. Magnetic materials 

Next we consider the case of magnetic materials, i.e. u r ^ 
1. In general, these materials are lossy and u' has relatively 
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Fig. 20. The simulated spectra of the complex permittivity and permeability 
for a lossy dielectric material with e r = 5 — 0.2i, fi r = 1. The upper graph 
presents the errors due to S-parameters, with ^SqI = —0.001 and 86 a = 
— 1°, the lower graph is simulated with 8L a i r = 0.1 mm. 

low values. 

The magnetic energy loss is included by assuming a non- 
zero value for the imaginary part of magnetic permeability. 
The simulations are done for three different materials with [i r 
= 2 — O.li, 2 — 0.3i and 2 — 0.5i. The dielectric loss is taken 
zero with e r = 5 + 0i. Fig. l2"Tlshows how the measured results 
would look like when only the L a i r error is considered. These 
plots were made using both n and z, e r = n/z and /i r = nz. 
In principle, it is possible to extract /i r as /i r = n 2 /e r , which 
would lead, much like in the case of dielectric materials, to 
lower errors. However, the typical situation in practice is that 
the experimentalist has no reliable knowledge of either e r or 
\i r for such materials. Similar to the lossy dielectric material 
previously simulated, the calculations show that there is no 
divergence in the permittivity and permeability spectra of lossy 
magnetic materials. 

VI. Conclusions 

This paper shows that the electromagnetic properties of 
materials, such as the complex refractive index, the complex 
permittivity and the complex permeability, can be measured 
using the transmission/reflection method without having to 
know the positions of the reference planes. The method is 
explicit (i.e. non-iterative) and easy to implement. The value 
of magnetic permeability is not assumed: both magnetic and 
non-magnetic material can be characterized. 
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Fig. 21. Simulated measurements of complex electromagnetic parameters for 
lossy magnetic materials with t r = 5, fi r = 2 — O.li (solid-blue), 2 — 0.3i 
(dashed-green) and 2 — 0.5i (dot-dashed-red). The error due to L a i r is set 
to 0.1 mm. 



The algorithm has been verified experimentally on two types 
of low-loss dielectric materials. The results show that at Fabry- 



Perot resonances, discontinuities occur in the real-part spec- 
tra of constitutive parameters. However, these discontinuities 
occur only over a few measured points, therefore in many 
situations of practical interest the results at these points can 
be neglected and the true values of material parameters can 
be found by interpolation. For the spectra of imaginary values, 
the method leads to results which are relatively stable all over 
the measurement bandwidth. 

The numerical analysis and the experiments show that the 
method has relatively low uncertainties at higher frequencies. 
This characteristic is an intrinsic property of the transmis- 
sion/reflection method because a significant change in the 
phase of the S-parameter due to the presence of the material 
requires that the sample length should be above a certain size, 
compared to the wavelength of the signal. 
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